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Abstract. This talk reviews the status of QCD calculations of photon and dilepton production rates 
in a Quark-Gluon plasma. Theses rates are known to order (D(a a ). Their calculations involve various 
resummations to account for well identified physical effects that are briefly described. Lattice calculations 
of the spectral functions give also access to the dilepton rates. Comparison with perturbative results 
points to inconsistencies in both approaches when the dilepton energy becomes small. 
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1 Introduction 

Photons or lepton pairs are produced at various stages 
of a nucleus- nucleus collision. Prompt photons and large 
mass dileptons are produced in the initial partonic col- 
lisions. Their rates can be calculated using zero temper- 
ature perturbative QCD. They populate the high energy 
part of the spectrum. All other photons or dileptons result 
from secondary interactions between the produced parti- 
cles. We focus here on the photons which are produced in 
a thermalized quark-gluon plasma (QGP). Their rates can 
be calculated using equilibrium thermal held theory. We 
shall not discuss how the rates can be combined with the 
space time evolution of nucleus- nucleus collisions in order 
to obtain the observed yields jUGlCIE])- Nor shall we dis- 
cuss photons produced in the hadronic phase (see the talk 
by Charles Gale in these proceedings). This talk builds 
on Ref. and extends some of the discussions presented 
there. 

The photon production rate can be expressed in terms 
of the the current-current correlator (j fJ _(0)j l ,(x)) , where 
the electromagnetic current is j^(x) — ip( x )^ffj.' l P( x )- To 
leading order in the electromagnetic fine structure con- 
stant a, the photon production rate reads [HIE] (uj 2 = q 2 )- 
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where n< u (tu,q) is the electromagnetic polarization ten- 
sor: 



(2) 



The second of eqs. Q gives the photon production rate 
in terms of the retarded polarization tensor. A similar for- 
mula exists for lepton pairs (Q = (ui, q), Q 2 = uj 2 — q 2 > 
0): 
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indicates a threshold at Q 2 = 4m 2 , with mi the mass of 
the lepton. 

In the first part of this talk, we shall review the an- 
alytical calculations of the rate, based on weak coupling 
techniques. Then we shall briefly discuss the estimates ob- 
tained from lattice determinations of spectral functions. 



2 Weak coupling calculations 
2.1 Leading order 



Fig. 1. Real processes contributing to photon and dilepton 
production up to 0(ot s ). One loop virtual corrections to the 
first process also contribute at this order. 
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The leading order contribution to the dilepton rate 
is obtained from the one-loop contribution to the polar- 
ization tensor, and corresponds to the Drell-Yan process 
illustrated by the diagram on the left of fig. ^ (only the 
production of the virtual photon is represented). It was 
evaluated for a QGP in 



2.2 First perturbative corrections 

The corrections of order 0(a s ), where a s = g 2 /4ir, with 
g the QCD gauge coupling, correspond to the two dia- 
grams in the right of figure ^ Their calculation reveals 
two problems. For the dilepton rate (Q 2 > 0) in a plasma 
of massless quarks and gluons, each individual contribu- 
tion to eq. J3J contains a mass singularity, and it is only 
after a careful summation of all the real and virtual cor- 
rections that one gets a finite result |1()U11II12| . This is 
nothing but a manifestation of the KLN theorem |14II13) . 
In the case of real photons (Q 2 — ► + ) a new singularity 
appears, with contributions of the form 



Im n Tet {u) , q) oc aa B \a(uT/Q 2 ) 



(5) 



at small Q 2 . The singularity originates from the presence 
of intermediate massless quarks (the "vertical" propaga- 
tors in the two diagrams in the right of fig-CJ. As we shall 
see, plasma effects induce effective masses and cure part 
of the difficulty. 



2.3 Scales and degrees of freedom in a quark-gluon 
plasma 

At this point it is useful to recall some basic properties of 
a quark-gluon plasma in the weak coupling regime, i.e. at 
sufficiently high temperature |15|. This regime is charac- 
terized by a hierarchy of momentum scales. Most of the 
plasma particles have momenta of the order of the tem- 
perature T, and since their density is of order T 3 , T is also 
the inverse of the inter-particle distance. Besides, collec- 
tive excitations can develop in the system. Such collective 
phenomena are particularly important at the scale gT, 
where g is the gauge coupling (the reason why these ex- 
citations are called collective is that, when j « 1, their 
wavelength ~ 1 / gT is large compared to the inter-particle 
distance ~ l/T, so that many particles participate in the 
excitation). Systematic corrections to the propagation and 
interactions of such collective excitations involve the re- 
summation of the so-called "hard thermal loops" |lf>U17| . 
Note that the soft collective modes also modify the spec- 
trum of hard particles, giving them a mass that we shall 
refer to ctS 777 oo (~ gT). Finally, another scale plays an im- 
portant role in a quark-gluon plasma: this is the scale g 2 T 
where perturbation theory breaks down because of the 
presence of unscreened magnetic fluctuations. The scale 
g 2 T characterizes also the rate of collisions with small 
(~ gT) momentum transfer. To see that write the scat- 
tering cross section as a = J d<7 2 (der/d<7 2 ), where typically 
da/dq 2 ~ <? 4 /? 4 - The collision rate is 7 = na, so that, 



with n ~ T 3 , 7 ~ <? 4 T 3 J dq 2 / q . The infrared divergence 
of the integral is cut-off by the screening mass mo (~ gT) 
(mo is an example of a "hard thermal loop") leaving a 
finite result 7 ~ g 4 T 3 /m D ~ g 2 T. This simple estimate 
applies when collisions involve dominantly small momen- 
tum transfer ~ gT. However, when calculating the effect 
of collisions on transport properties, the dominant colli- 
sions involve large angle scattering and the infrared cut-off 
is actually taking place at a larger momentum scale, of or- 
der T. Thus, most transport coefficients end up being of 
order g A T\n(l/g) [Hj. 



2.4 Resummation of hard thermal loops 

We can now return to eq. © and observe that the log- 
arithmic singularity at Q 2 — > is due to the exchange 
of a soft massless quark. Once the HTL correction is in- 
cluded on the quark propagator, the quark effectively ac- 
quires a mass moo of order gT ( m 2 ^ = ira s CfT 2 with 
Cf = (N 2 — 1)/2N C ). By taking into account this thermal 
correction one obtains a finite photon polarization tensor 
|19U2()| . For hard photons, it reads: 



Im77 re t / V(w,<?) 
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The numerical factor 5/9 is the sum of the quark electric 
charges squared for 2 flavors (u and d); for 3 flavors (u, d 
and s), this factor should be replaced by 6/9. This formula 
indicates how the infrared problem is cured: Q 2 in Eq. (JjjJ 
is effectively replaced by m 2 ^ in the logarithm as soon as 
Q 2 becomes small compared to m 2 ^. 




< 



Fig. 2. Processes that are promoted to 0(a s 
singularities. 



by collinear 



This, however, is not the final answer for the photon 
and dilepton rates at 0(a s ). Indeed there are other pro- 
cesses, formally of higher order, which are strongly en- 
hanced by collinear singularities and become effectively of 
order a s . This was first realized for soft photon production 
by quark bremsstrahlung 21,22 (third diagram of fig. [21 
starting from the left). The diagram on the right of figure 
121 shares the same property, but contributes significantly 
only to hard photon production [22], due to phase-space 
suppression in the case of soft photons. Note that a naive 
power counting would indicate that these two diagrams 
contribute to 0(a 2 ). These two diagrams represent colli- 
sion processes: they originate from cutting a loop insertion 
in the gluon propagator in order to get the imaginary part 
(the cuts are indicated by the dotted lines in the two dia- 
grams on the left of fig. El- 

In order to understand the origin of the "collinear en- 
hancement" , let us focus on the quark propagator between 
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P+Q 




Fig. 3. The virtual quark of momentum P + Q emitting a 
real photon (Q 2 — 0) and an on-shell quark of momentum P 
(P 2 = m 2 ). 



the quark-gluon vertex and the photon emission. The vir- 
tually of this off-shell quark is easily estimated (see fig-EJ : 



(P + Q) 2 



2P-Q 



Pz 



where rtv\ 



A 



m;^. Thus the virtuality of this quark 
can become very small if the quark is massless and the 
photon is emitted forward ( p± — > 0). Of course, the 
quark thermal mass m^ prevents these diagrams from 
being truly singular. However, contrary to the 0(a s ) dia- 
grams, the singularity is here linear instead of logarithmic, 
and brings a factor T 2 /m 2 ^ ~ l/a s . Combined with the 
a 2 that comes from the vertices, the singularity turns the 
contribution of these diagrams into an order 0(a s ) con- 
tribution. This was evaluated in |22] and [231 ( see also [21] 
and where an erroneous factor 4 was pointed out), and 
a closed expression was obtained in [2323- The result is 
of the form: 



Im/Jret^lw, q) = const aa s 



T 3 

IT 2 +UJT 

UJ 



(8) 



In this formula, the term in 1 jio dominates for soft photons 
and comes from the bremsstrahlung diagram, while the 
term in ui comes from the second diagram and dominates 
the rate of very hard photons (w » T). 

It is worth mentioning that the purely numerical pref- 
actor (not written explicitly in the previous formula) is 
a function of the ratio of the quark thermal mass 
to the gluon Debye mass m D . In the HTL approxima- 
tion, this ratio is a constant independent of the coupling 
and temperature, that depends only on the number of col- 
ors and flavors; for 3 colors and Nf flavors, this ratio is 
moc /m D = y/2/{6 + N f ). 

This enhancement due to a quasi-collinear emission of 
the photon also occurs for the emission of virtual photons 



with a small invariant mass (Q 



q 2 < uj 2 ), but be- 



comes less and less important when the photon invariant 
mass increases. For virtual photons with vanishing mo- 
mentum (i.e. for which the invariant mass is maximal, at 
a given energy) , the two-loop diagrams contribute only at 
the order g 3 [28]. instead of g 2 for real photons. 



2.5 LPM resummation 

The collinear enhancement that we have identified on so- 
me of the order a 2 processes affects in fact an infinite set of 
processes. In order to explain the issue in physical terms, it 
is convenient to introduce the concept of photon formation 



time. Let us return to the process in figure [31 The photon 
formation time can be identified with the lifetime of the 
virtual quark, which is itself related to its virtuality by 
the uncertainty principle. A simple calculation gives: 



SE 



q 



E„ — E, 



p+q 



2 p z (p z +uj) 



(9) 



where the 3-momentum of the photon defines the longitu- 
dinal axis. The formation time is t F = 1/SE. The collinear 
enhancement in the diagrams of figure [21 due to the small 
virtuality of the quark that emits the photon, can be 
rephrased in terms of the large photon formation time. 



(7) Typically, in a quark-gluon plasma, we have m\ ~ g 2 T 



while p z ~ T, so that SE in eq. © is SE ~ g 2 T for ui ~ T. 
That is, the photon formation time is of the same order 
of magnitude or larger than the quark mean free path be- 
tween two soft collisions, i.e. t F ~ I/7, where 7 ~ g 2 T is 
the collision rate estimated earlier. Note that the estimate 
done earlier is indeed the relevant collision time scale for 
the production of photons almost collinear with the charge 
particle, that is with a typical transverse momentum of 
order gT: this is the kinematical condition leading to the 
enhancement that we are discussing (by "enhancement", 
we mean, as earlier, the phenomenon by which higher or- 
der diagrams turn out to contribute at the same order as 
a given elementary process). The sensitivity to the colli- 
sional width found in 2H| occurs in the same kinematical 
conditions. When the mean free path becomes of the order 
of, or smaller than the photon formation time, the effects 
of multiple collisions on the production process can no 
longer be ignored. The result of such multiple scattering 
is to reduce the rate compared to what it would be if all 
collisions could be treated as independent source of pho- 
ton production. This phenomenon is known as the Landau 
Pomeranchuk Migdal (LPM) effect GHEES]- 



1" 



< 



Fig. 4. Diagram illustrating the conditions under which mul- 
tiple collisions need to be taken into account in the photon or 
dilepton production process: / is the typical range of the inter- 
action producing the collisions; A is the mean free path, and 
tF the photon formation time. 



While the early treatment of the multiple scattering 
was done in terms of kinetic equations, modern discussions 
used the language of quantum field theory. The multiple 
scattering diagrams that must be resummed are the lad- 
der diagrams, as illustrated in figure [3] (these are the typ- 
ical diagrams that are taken into account by a Boltzmann 
equation |33|L Cancellations between self-energy correc- 
tions and vertex corrections remove any sensitivity to the 
magnetic scale: physically, such cancellations reflect the 
fact that ultrasoft scatterings (at momentum scale softer 
than gT) are not efficient enough to produce a photon. A 
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thorough diagrammatic analysis explaining why it is the 
ladder family of diagrams that needs to be resummed in 
order to obtain the complete leading 0(a s ) photon rate is 
presented in |34| . 



2.6 Some numerical results 

The integral equation was solved numerically in |35|. and 
the results are displayed in figure In this plot, 'LPM' 




a =0.3, 3 colors, 3 flavors, T=l GeV 



Fig. 5. Rcsummation of ladder diagrams. 

In the recent literature, the resummation of the ladder 
diagrams is presented as follows. The photon polarization 
tensor is written explicitly as |34U35U3b| : 



dpo [n F {r ) - n F (p )] 

-OO 

d 2 P± 



pi 



x Re 



P± ' fiP±) , 



iporo) 2 
(10) 



with r = po +qj, n F (p ) = 1 / (exp (p /T) + 1) the Fermi- 
Dirac statistical weight, and where the dimensionless func- 
tion /(p_|_) denotes the resummed vertex connecting the 
quark line and the transverse modes of the photon 1 . This 
function is dotted into a bare vertex, which is proportional 
to p ± . The equation that determines the value of f(p±) 
is a Bethe-Salpeter equation that resums all the ladder 
corrections to the vertex 34 .3311361 : 



-/(Pj 



2Pj 



4ira s C f T 



d 2 l ± 



■C(Ix)[/(Px + I±) 



f(p±)) , 
(11) 

where t F is the formation time, t F = 1/SE, with SE given 
by eq. © , and where the collision kernel has the following 
expression |2fi| : 

C(ix) = A- 7X J -^- (12) 



where the two terms correspond to the exchange of a trans- 
verse and a longitudinal gluon, respectively. Note that the 
quark propagators should be dressed in a way compatible 
with the resummation performed for the vertex, in order 
to preserve the gauge invariance: this is the origin of the 
term —f(p±) under the integral in eq. which has 

the effect of resumming the collisional width on the quark 
propagator. From this integral equation, it is easy to see 
that each extra rung in the ladder contributes a correction 
of order a s Tt F ~ 0(1) since t F ~ l/g 2 T. Therefore, all 
these corrections contribute to 0(a s ) to the photon rate. 
Note again that the only parameters of the QGP that en- 
ter this equation are the quark thermal mass and the 
Debye screening mass m D . 

For the emission of real (massless) photons, only the trans- 
verse polarizations of the photon matter. 
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Fig. 6. 0(a s ) contributions to the photon production rate 
in a QGP. The parameters used in this plot are a a = 0.3, 3 
colors, 2 flavors and T = 1 GeV. From 35 . 



denotes the contribution of all the multiple scattering di- 
agrams, while '2 — ► 2' denotes the processes of figure ^ 
The single scattering diagrams (figure [5J are also given so 
that one can appreciate the suppression due to the LPM 
effect (ranging typically from 15 to 30%). 

Dilepton production basically suffers from the same 
problems, and the solution follows the same path. Two dif- 
ferences are worth mentioning here. First of all, the Drcll- 
Yan process qq — > 7* — > contributes if Q 2 > Aro^. 

The Drell-Yan process has been evaluated in [H], the 2-^2 
processes have been evaluated in |12j . In addition, vir- 
tual photons have a physical longitudinal mode that con- 
tributes to the rate of lepton pairs. In order to take this 
mode into account, one must introduce a scalar function 
g(p±) similar to f(p±_), which describes the coupling of 
the quark line to a longitudinal photon. This new vertex 
function obeys an integral equation |37| similar to eq. Ijlll) . 
that resums the corrections due to multiple scatterings. 
This new integral equation can also be solved numerically, 
and the resulting dilepton rate (for the same parameters as 
in figureEland a total energy of the pair set to ui = 5 GeV) 
is plotted in figured One can see that the multiple scat- 
tering corrections are important for all pair masses below 
the threshold of the Drell-Yan process. Note also that the 
threshold of the tree-level process is completely washed 
out when multiple rescatterings are resummed. 



3 Lattice calculations 

Attempts to calculate directly on the lattice the produc- 
tion rate of dileptons in a quark-gluon plasma appeared a 
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Fig. 7. 0(a s ) contributions to the dilepton production rate Fig. 8. Lattice result for the production rate of dileptons with 



in a QGP. From |SZ|. 



q) = 0. From [H]. 



few years ago [33] • What can be calculated on the lat- 
tice is the Euclidean correlator of two vector currents, 
7J(r, x) = (j M (0, 0)j^(r, x)), where r G [0, 1/T] is the Eu- 
clidean time. It is also easy to obtain the spatial Fourier 
transform at zero momentum, 77(r, q — 0), by just sum- 
ming over the spatial lattice sites. The imaginary part of 
the real time self-energy is then related to this object by 
a simple spectral representation: 



n{r,q) 



°dto 

ImiZret /J 



cosh(w(r- 1/2T)) 
sinh(w/2T) 



(13) 

This equation uniquely defines Im/J ret %(a;, q) if il(r, q) 
is known for all r G [0, and if one prescribes the 
behavior of the solution at large to. 

However, the function 77(r, q = 0) is known only on 
the discrete temporal lattice sites, which prevents us from 
determining uniquely Im7T rc t M ^(w, 0). This problem has 
been reconsidered recently using the Maximum Entropy 
Method |38II39| , which is a way to take into account prior 
knowledge about the solution (positivity, behavior at the 
origin, etc..) in order to determine the most probable so- 
lution compatible with the lattice data and with this a 
priori information. The result obtained for zero momen- 
tum dileptons via this method is displayed in Fig. [H] for 
two different values of the temperature. Note that this is 
a quenched lattice simulation. This result displays several 
interesting properties. At energies above 4T, the full rate 
is very close to the contribution of the Born term, while at 
energies smaller than 3T it drops to extremely small val- 
ues. In addition, when plotted against to/T, the curves for 
the two temperatures fall almost on top of one another, 
indicating that the result scales like a universal function 
of to/T, at least within the errors. 

The suppression at small to has attracted a lot of inter- 
est because it contradicts expectations based on pertur- 
bation theory: the resummation of thermal masses would 
indeed produce a drop of the Born term because of thresh- 
old effects (sec Eq. (QJ), but higher order processes that 
do not have a threshold would fill the spectrum at small 



to. Also, a threshold related to the quark masses would 
occur at much smaller to than 3T, since thermal masses 
are typically m ~ gT. Finally one may question whether 
the accuracy of lattice calculations in the small to regime 
is not spoiled by finite volume effects. 

On the other hand the polarization tensor at small 
frequency is related to the electric conductivity 0U| by 
the relation: 



(T e i = lim Im/J TO ti( w ;0)/ 6w - 

uj— >0 



(14) 



From this relation, one expects Im7T rct *(w, q = 0) cx to 
when to — > 0. This implies that the static dilepton rate 
should diverge when to — ► 0. Unless the electric conduc- 
tivity in quenched QCD is nearly zero for some reason, 
the lattice dilepton rate disagrees with this prediction at 
small to. Note that 'small' in these considerations means 
an to small enough to be in the hydrodynamical regime, 
i.e. to < g 4 T. In a strong coupling theory, this regime could 
start as early as to ~ T. 

Note that the previous argument rests on the possi- 
bility to replace Im77 rct Al A1 (w, 0) (involved in the calcu- 
lation of the dilepton rate) by ImII Te t'i(to,0) (which is 
the quantity needed to calculate the conductivity). This 
is guaranteed by the Ward identity 7J ret A1 A1 (w, 0) = 0. 
From this it follows indeed that, unless singularities occur, 
77 rot °V,0) = 0. 

The electric conductivity has been calculated on the 
lattice by S. Gupta 0O|> and a finite result was obtained. 
This calculation provides an illustration of the sensitivity 
of the maximum entropy method used to reconstruct the 
spectral functions to the prior information. Gupta's calcu- 
lation assumes explicitly that the spectral function that he 
wants to determine behaves linearly in to at small to. The 
maximum entropy procedure yields then a finite value for 
the slope, i.e., a finite value for the electric conductivity. In 
|38| on the other hand, no particular assumption is made 
about the behavior of the spectral function at small to. For 
further discussion on the difficulty of extracting transport 
coefficients from Euclidean lattice correlators, see |4 11142] 
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4 Discussion and outlook 

As of now, there are in fact arguments indicating that both 
the perturbative calculations and the lattice calculation 
are incorrect at small ui. If one evaluates eq. at r = 
1/2T, one gets a sum rule: 



Im77 rct %(co,g = 0) 
sinh(w/2T) 



= 77(1/27; 0) < oo 



(15) 



This sum rule is violated by all the existing analytical 
weak coupling calculations (they give an infinite result). 
For instance, the expression in eq. (jHJl of the imaginary 
part of 77 diverges as 1/lo at small u>. The LPM effect 
would reduce the divergence to one in 1/ y/U. But the ex- 
pected linear behavior is not achieved in present approxi- 
mations. In fact, none of the existing calculations includes 
correctly the dissipative effects that appear when one en- 
ters the hydrodynamical regime (uj — *■ 0). 

As for the lattice, consider again eq. (|15fl. If one as- 
sumes that the integral is dominated by the behavior of 
Im 77 at small to, i.e. Im 77* (a;, q = 0) ~ 6cr c i^, one obtains 
an estimate of the integral 



7Zj(r = l/2T,g = 0)« / -- 



sinh(w/2T) 



= 6<j c \itT 



(16) 



One could perhaps argue that this estimate provides in 
fact an upper bound for 77* (r = 1/2T, q — 0) (if we admit 
that the actual function never exceeds the linear extrap- 
olation). In order to make contact with lattice estimates, 



we define C EM = 
equation above as 



Aira^r eh. Then we can rewrite the 



1 77] (r = 1/2T,0) 



(17) 



The left hand side can be obtained from Gupta's calcu- 
lation and is a number of order 7. The right hand side 
is given in and is weakly dependent of the tempera- 
ture; it is a number of order 0.12. These simple estimates 
suggest that the lattice calculations are not fully consis- 
tent. If we would admit that eq. I|16(l provides a lower 
bound for the electric conductivity, then the calculation 
in |38j should yield a finite a e i, which is not compatible 
with fig. |H1 It is also somewhat puzzling that the value of 
a obtained in 00] is so much larger than the simple esti- 
mate based on eq. (|16|) : it would be interesting to know 
whether the values of 77(1/2T,0) obtained in gO] agree 
with those given in |ri8| . 
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